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ABSTRACT 

We discuss the origin of the very different hard X-ray spectral shapes observed in the 
Low Hard State (LHS) and High Soft State (HSS) of accreting black holes. We study 
the effects of synchrotron self- absorption on the Comptonising electron distribution in 
the magnetised corona of accreting black holes. We solve the kinetic equations assum- 
ing that power is supplied to the coronal electrons through Coulomb collisions with a 
population of hot protons and/or through the injection of non-thermal energetic elec- 
trons by some unspecified acceleration process. We compute numerically the steady 
state particle distributions and escaping photon spectra. These numerical simulations 
confirm that synchrotron self-absorption, together with e-e Coulomb collisions consti- 
tute an efficient thermalising mechanism for the coronal electrons. When compared 
to the data they allow us to constrain the magnetic field and temperature of the hot 
protons in the corona independently of any dynamical accretion flow model or ge- 
ometry. A preliminary comparison with the LHS spectrum of Cygnus X-1 indicates a 
magnetic field below equipartition with radiation, suggesting that the corona is not 
powered through magnetic field dissipation (as assumed in most accretion disc corona 
models). However, in the LHS of Cygnus X-1 and other sources, our results also point 
toward proton temperatures lower than 3 x 10^*^ K, i.e. substantially lower than typ- 
ical temperatures predicted by the ADAF-like models. In contrast, in the HSS both 
the proton temperature and magnetic field could be much higher. We also show that 
in both spectral states the magnetised corona could be powered essentially through 
acceleration of non-thermal particles. Therefore, contrary to current beliefs, energy 
dissipation processes in the corona of the HSS and that of the LHS could be of very 
similar nature. The main differences between the LHS and HSS coronal emission can 
then be understood as the consequence of the much stronger radiative cooling in the 
HSS caused by the soft thermal radiation coming from the geometrically thin accretion 
disc. In the LHS the soft cooling photon flux is much weaker because the accretion 
disc is either truncated at large distances from the black hole, or much colder than in 
the HSS. 

Key words: accretion, accretion discs - black hole physics - radiation mechanisms: 
non-thermal - methods: numerical - gamma-rays: theory - X-rays: binaries 



1 INTRODUCTION 

Black hole binaries are observed in two main spectral states. 
At luminosities exceeding a few percent of the Eddington lu- 
minositjQ (Le), the spectrum is generally dominated by a 
thermal component peaking at a few keV which is believed 



E-mail:malzac@cesr.fr 
^ For hydrogen gas, the Eddington luminosity is Le = 1.3 x 
10^^ M erg s~^ where M is the black hole mass expressed in units 
of solar masses 



to be the signature of a geometrically thin optically thick ac- 
cretion disc (Shakura & Sunyaev 1973). At higher energies 
the spectrum is non-thermal and usually presents a weak 
steep power-law component (photon index F ~ 2.3 — 3.5) 
extending at least to MeV energies, without any hint for 
a high energy cut-off. This soft power law is generally in- 
terpreted as inverse Compton up-scattering of soft photons 
(UV, soft X) by an hybrid thermal/non-thermal distribu- 
tion of electrons in a hot relativistic plasma (the so-called 
"corona"). The electron distribution is then characterised 
by a temperature ~ 30-50 keV, a Thomson optical depth 
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~ 0.1-0.3 and a quasi power-law tail (71(7) oc 7~^, with in- 
dex s in the range 3.5-4.5) which is responsible for most of 
the hard X-ray and 7-ray emission (Gierlihski et al. 1999, 
hereafter G99; Frontera et al. 2001a; McConnell et al. 2002, 
hereafter MC02; Del Santo et al. 2008). Since in this state 
the source is bright in soft X-rays and soft in hard X-rays it 
is called the High Soft State (hereafter HSS). 

At lower luminosities (L< 0.01 Lb), the appearance of 
the accretion flow is very different: the X-ray spectrum can 
be modelled as a hard power-law F ~ 1.5 — 1.9 with a cut- 
off at ~ 100 keV. The vFu spectrum then peaks at around 
a hundred keV. Since the soft X-ray luminosity is faint and 
the spectrum is hard, this state is called the Low Hard State 
(hereafter LHS). LHS spectra are generally very well fit- 
ted by thermal Comptonisation, i.e. multiple Compton up- 
scattering of soft photons by a Maxwellian distribution of 
electrons (see Sunyaev & Titarchuk 1980) in a hot (fcTc ~ 
50-100 keV) plasma of Thomson optical depth tt of order 
unity. Although a weak non-thermal component is clearly 
detected at MeV energies in the LHS of Cygnus X-1 (MC02) 
and GX 339-4 (Wardzihski et al. 2002; Joinet et al. 2007) 
most of the hard X-ray luminosity emerges in the form of 
Comptonisation by a thermal electron distribution. Spectral 
fits with hybrid thermal/non-thermal Comptonisation mod- 
els suggest the slope of the non-thermal tail in the electron 
distribution is steeper (s > 5) than in the HSS (see MC02). 

In brief, both spectral states are well represented by 
Comptonisation by an hybrid electron distribution. In the 
LHS the temperature and optical depth of the thermal elec- 
trons are higher, and the slope of the non-thermal tail seem 
steeper than in the HSS. Consequently, the X-ray emission 
is dominated by thermal Comptonisation in the LHS and by 
non-thermal Comptonisation in the HSS. 

The different spectral states are usually understood in 
terms of changes in the geometry of the accretion fiow. Ac- 
cording to a popular scenario (see e.g. Done, Gierlihski and 
Kubota 2007, hereafter DGK07), in the HSS a geometrically 
thin accretion disc extends down to the last stable orbit and 
is responsible for the dominant thermal emission. This disc is 
the source of soft seed photons for Comptonisation in small 
active coronal regions located above and below the disc. The 
magnetic field lines rise above the accretion disc through 
magnetic buoyancy, transporting a significant fraction of the 
accretion power into the corona where it is then dissipated 
through magnetic reconnection (Galeev, Rosner, & Vaiana, 
1979). This leads to particle acceleration in the corona. A 
population of high energy electrons is formed which then 
cools down by up-scattering the soft photons emitted by the 
disc. This produces the high energy non-thermal emission 
(see e.g. G99; Zdziarski et al. 2002). 

In the LHS, the standard geometrically thin disc does 
not extend to the last stable orbit. Instead, the weakness 
of the thermal features suggests that it is truncated at dis- 
tances ranging from a few tens to a few thousand gravita- 
tional radii from the black hole (typically 1000-10000 km). 
In its inner parts, the accretion flow takes the form of a hot 
geometrically thick, optically thin disc (Esin, McClintock 
and Narayan 1997; Poutanen, Krolik and Ryde 1997). In 
the standard hot flow solutions (see e.g. Shapiro, Lightman 
and Eardley, 1976; Ichimaru 1977; Narayan and Yi 1994; 
Abramowicz et al 1996, Blandford and Begelman 1999), the 



gravitational energy is converted through the process of vis- 
cous dissipation into the thermal energy of ions. 

The ions only cool by heating up electrons through 
Coulomb collisions. This process being rather inefficient it 
leads typically to high proton temperatures (Ti ~ 10^^ K). 
Such electron heating naturally forms nearly Maxwellian dis- 
tributions of electrons (Nayakshin & Melia 1998). These 
electrons can then radiate their energy by Comptonising 
the soft photons coming from the external geometrically 
thin disc, as well as IR-optical photons internally generated 
through self-absorbed synchrotron radiation. The balance 
between heating and cooling determines the electron tem- 
perature which turns out to be close to 10^ K, as required to 
fit the LHS spectra (but see Yuan & Zdziarski 2004 and Sec- 
tion [5]3}. These two-temperature accretion flows are domi- 
nated by advection over a wide range of the parameter space, 
they are usually referred to as ADAF (Advection Dominated 
Accretion Flows). 

However, Coulomb heating is the dominant electron 
heating process only if the magnetic fleld is strongly sub- 
equipartition (ratios of gas to magnetic pressure > 10). For 
larger magnetic fields, turbulent stresses may heat the elec- 
trons directly in the same way as they heat the protons. For 
approximately equipartition magnetic field the turbulence 
primarily heats the electrons (Quataert and Gruzinov 1999). 
For even larger magnetic fields reconnection may be the 
dominant channel of energy dissipation into electron heat- 
ing (Binovatyi-Kogan and Lovelace 1997). From an obser- 
vational perspective, recent spectral modelling of the hard 
X-ray spectra of Sagitarius A* and black holes binaries in 
the LHS suggest that a significant fraction (~50 percent) 
of the energy is dissipated directly through electron heating 
(Yuan et al. 2003; Yuan et al. 2005; 2007) possibly through 
MHD turbulence, reconnection, and weak shocks. Such pro- 
cesses are likely to generate a non-thermal electron distribu- 
tion function at high energy (see e.g. Li & Miller 1997). If 
this is the case and if heating processes are essentially non- 
thermal in both spectral states, we have to understand the 
cause of the observed difference in the distribution of the 
Comptonising electrons (i.e. thermal LHS and non-thermal 
in the HSS). 

The main alternative to the ADAF-like models for the 
LHS relies on electron heating by magnetic reconnection 
above the standard thin disk (Liang & Price 1977; Galeev 
et al 1979) as discussed earlier for the HSS. Since the pro- 
cess of magnetic reconnection remains poorly understood, 
this model has not yet been elaborated in detail, although 
numerical simulations have shed some light on this problem 
(Miller & Stone 2000; Hirose, Krolik & Stone 2006). However 
the free parameters of this accretion disc corona model (see 
Haardt, Maraschi & Ghisellini 1993) are very constrained 
by the observations (see Malzac 2007 for a discussion). This 
model successfully reproduces LHS spectra by postulating 
the presence of outfiowing magnetic active region above a 
cold thin disc extending close to the last stable orbit (Be- 
loborodov 1999; Malzac, Beloborodov & Poutanen 2001). 
In this framework, the thin disc is never truncated, simply 
colder in the LHS (see Malzac 2007). This dynamic corona 
is assumed to be powered through the same magnetic pro- 
cesses as the accretion disc corona of the HSS. However the 
cause for the different electron distributions was never dis- 
cussed so far, which is precisely the aim of this paper. 
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The main idea we investigate here is that in presence 
of a magnetic field, the electron distribution can appear 
thermal even when acceleration mechanisms would produce 
non thermal distributions. Synchrotron self-absorption in- 
deed represents a very fast and efficient thermalising mech- 
anism. The importance of this effect for the emission of 
compact sources was first pointed out by Ghisellini, Guil- 
bert and Svensson (1988). These authors and then Ghis- 
ellini, Haardt and Svensson (1998) considered the case of 
high energy electrons continuously injected in a magnetised 
region and then cooling down through synchrotron and in- 
verse Compton emission. Solving the kinetic equations, they 
showed that due to the very fast emission and absorption 
of synchrotron photons, the electron distribution relaxes to- 
ward a Majcwellian distribution in a few light crossing times. 
They named this effect the synchrotron boiler. The steady 
state distribution then consists in a Maxwellian with a quasi- 
powerlaw high energy tail (formed by the cooling electrons) . 
In the context of ADAF models, Mahadevan and Quataert 

(1997) showed that the synchrotron boiler together with e-e 
Coulomb collisions would efficiently thermalise the electron 
distribution of the hot accretion flow. 

Although these studies demonstrated the correct quali- 
tative behaviour, they were based on a very simplified treat- 
ment, in which the electron and photon kinetic equations are 
not solved self-consistently and the Compton process is mod- 
elled in the Thomson approximation. Although this treat- 
ment provides important insights, it is not realistic enough 
for a quantitative comparison with data. Recently, Vurm and 
Poutanen (2008), presented the results of a computer code 
which overcomes these problems and confirms the earlier 
results. Their code solves simultaneously the kinetic equa- 
tions for photons and electrons using the full Klein-Nishima 
cross section for Compton scattering. Their simulations how- 
ever do not account for other radiation processes such as 
pair production, or other, possibly important, thermalising 
mechanisms such as Coulomb interactions between elec- 
trons. Several other computer codes such as eqpair (Coppi 
1992, 1999) or the Fokker-Planck code by Nayakshin & Melia 

(1998) take into account these processes, but not the ther- 
malising effects of the magnetic field. Moreover, eqpair pos- 
tulates a priori a low energy thermal electron distribution, 
which is not ideal when studying thermalisation processes. 

The first full treatment of the problem (including all the 
relevant processes) under the 1-zone, homogeneous, isotropic 
approximation, was performed only recently with a code pre- 
sented in Belmont, Malzac & Marcowith (2008a, 2008b). In 
the present paper we use this code to investigate in more de- 
tails the synchrotron boiler effects in physical situations rel- 
evant to the corona of accreting black holes. We also present 
a first direct comparison of our simulations with the data. 

In Section [2] we present the assumptions of our model 
and derive useful analytical estimates for the temperature 
of the hot protons and the equi-partition magnetic field as a 
function of compactness. In Section[3]we explore the param- 
eter space. Assuming that electrons are accelerated accord- 
ing to a power-law distribution, we show that, for reasonable 
parameters, the steady state spectrum is very similar to that 
observed in the LHS of black hole binaries. Then we show 
that the change in the thermal vs non-thermal Compton 
emission between the LHS and IfSS can be understood as 
the effect of the strong soft photons fiux cooling down the 



thermal part of the electron distribution to temperatures 
at which its Compton emissivity is negligible and leaving a 
bare non-thermal spectrum. We also discuss the effects of 
Coulomb heating by hot protons. Finally, in Section |4] we 
attempt to model observed LHS and HSS spectra of Cygnus 
X-1. This allows us to put some constraints on the mag- 
netic field and proton temperature of the coronal plasma in 
both states. This comparison also confirms that both spec- 
tral states can be modelled with power-law electron acceler- 
ation as unique heating mechanism. These results are then 
summarised and discussed in Section [5] 



2 MODEL 

2.1 Framework 

We use the code of Belmont et al. (2008b). This code solves 
the kinetic equations for photons, electrons and positrons 
accounting for Compton scattering (using the full Klein- 
Nishima cross section), electron-positron pair production 
and annihilation. Coulomb interactions (electron-electron 
and electron-proton), synchrotron emission and absorption 
and e-p bremsstrahlung. Radiative transfer is dealt using a 
usual escape probability formalism. This code was critically 
tested and found consistent with the results of the eqpair 
code (Coppi 1992) and the Fokker-Planck code of Nayak- 
shin and Melia (1998) which however do not account for 
synchrotron. 

Our code is fully time-dependent. An example of time- 
dependent simulation of the synchrotron boiler effect is pre- 
sented in Belmont et al. (2008b). For compactness param- 
eters of order of unity, or larger, the system relaxes to- 
ward equilibrium very quickly (typically a few light-crossing 
times). This relaxation time is smaller than, (or at most 
comparable to) the dynamical time-scale of the accretion 
disc at all radijf). Moreover the observation shows that the 
X-ray variability of black hole binaries and active galactic 
nuclei is dominated by fluctuations on time-scales that are 
much longer than the light-crossing time of the inner disc 
region which is supposed to produce those X-rays (see e.g. 
DGK07). Therefore, for most applications assuming steady 
state should be a good approximation. In this paper we limit 
our study to steady-state particle distributions and photon 
spectra. 

We consider a sphere with radius R of fully ionised 
proton-electron magnetised plasma in steady state. The 
Thomson optical depth of the sphere is tt = n + Ts, where 
Ti — TiifTT-R is the optical depth of ionisation electrons (asso- 
ciated with protons of density rii) and Ts — 2n^+OTR is the 
optical depth of electrons and positrons due to pair produc- 
tion {n^+ is the positron number density), ot is the Thom- 
son cross section. The radiated power is quantified through 
the usual compactness parameter: 

j_ Lot _ L lO^cm _ L SORg 

Rm,c3 1037erg R IQ-^Le R ' ^ ' 

^ The light crossing time tiight of a-n emitting region of size R, is 
always much shorter than the local dynamical time-scale tdyn at 
a distance R from a black hole. For a black hole of gravitational 
radius Rq, tdyn/*iig ht = 2-k^R/Rg 
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where L is the luminosity of the Comptonising cloud, rric the 
electron rest mass and c the speed of light. We introduced 
the Eddington luminosity Le = 1.3 x 10* M/M© erg s"\ 
and the gravitational radius Rq = GM/c?, where G is the 
gravitational constant and M the mass of the black hole. 



2.2 Power supply 

We will consider three possible channels for the energy in- 
jection in our coupled electron-photon system: 

(i) Non-thermal electron acceleration with a compactness 
^nth- We model the acceleration process by assuming elec- 
trons are continuously injected with a power-law distribu- 
tion of index Finj (i.e. n(7) oc 7~'"'°-'), between Lorentz fac- 
tors 7min and 7max. Obscrvatious of Cygnus X-1 in the soft 
state indicate Finj in the range 2.0-3.5 (G99). In most of 
our simulations, we set 7min = 1, 7max = 1000. As will 
be shown in section I3TT1 as long as 7min < 2 and 7max ^ 1, 
and Finj > 2 the exact values have little effects on the result. 
In this paper, we assume that the particles are reaccelerated 
rather than injected, and therefore we balance the number of 
'injected' non-thermal particles by removing the same num- 
ber of particles from the steady state distribution. We define 
the non-thermal compactness as the net energy input due to 
the acceleration process, i.e: 



/nth = [{7)i - (7)]n<7T/(i?c) 



(2) 



where n is the total number of injected/accelerated particles 
per unit time, (7) is the average Lorentz factor of electrons 
in the hot plasma, and (7)1 is the average Lorentz factor 
at which they are accelerated/injected (which depends on 
7min, 7max and Finj). This definition of /nth corresponds 
to that of the version of eqpair used by MC02 and other 
authors to fit the spectra of black hole binaries. We note 
that this assumes that re-acceleration occurs uniformly at 
all energies up to 7max (while in some other acceleration 
models the accelerated particles are drawn from the thermal 
part of the distribution). 

(ii) Coulomb heating with a compactness Ic- Electrons are 
supposed to interact by Coulomb collisions with a distribu- 
tion of thermal protons. When the protons have a larger 
temperature than the electrons, the electrons gain energy 
and Zc > (see below Eq.[3|. In our model, we fix Ic, then de- 
termine the proton temperature corresponding to this com- 
pactness and to the exact steady state electron distribution. 

In the case of a pure Maxwellian plasma, the dependence 
of the Coulomb compactness on electron and proton tem- 
peratures can be approximated as follows (Dermer 1986): 



Zc = VSTrrTTi In A(6'i -OAle^ + OA 

mi 



-3/2 



(3) 



where In A ~ 20 is the Coulomb logarithm and 9c — is 
the electron temperature in units of the electron rest mass. 
Equation [3] is valid in the limit of small 6i and 8e- For typical 
values 0c ~ 0.18 and Oi ~ 0.01 , h ~ 12rTri. 

To the first order in 6i/6c equation Q can be rewritten 



1 + 7 



gic 

TTTi 



where g{9c) = V79c (— 1 in the LHS). The proton tem- 
peratures provided by equation Q agree with the values 
obtained from our numerical simulations within 30 percent, 
even when the lepton distribution strongly departs from 
a Maxwellian. The largest discrepancies occur when the 
approximation of small 9i/9c breaks down. For the typical 
electrons temperatures observed in accreting black hole 
sources, this happens for proton temperatures > 50 MeV. 
In sources with accurately measured electron temperature, 
Thomson optical depth, and luminosity, and for a given size 
of the emitting region, equation Q provides an upper limit 
on the proton temperature. It will be used in Sections |4] 
and 15.31 to constrain the temperature in Cyg X-1 and 
XTE J1118+480. 

(iii) External soft radiation coming from the geometri- 
cally thin accretion disc and entering the corona with a 
compactness Is, that we model as homogeneous injection of 
photons with a blackbody spectrum of temperature fcTbb. 

Since, in this model, all the injected power ends up into 
radiation we have: I = Znth + lc + L in steady state. 

We note that this set-up is very similar to the one used 
in the eqpair model when fitting the data of black hole bi- 
naries (see G99). The main difference is that we take into 
account the synchrotron emission and absorption processes 
and the associated synchrotron boiler effect. In our code 
the thermalisation process is treated self-consistently while 
in EQPAIR it is assumed that the electron distribution is 
Maxwellian at low energies. Another difference is that we 
have a detailed description for Coulomb heating by hot pro- 
tons. In EQPAIR, electron heating is parametrised through a 
'thermal' compactness Ith quantifying the amount of power 
provided as heating to the thermal pool of cool electrons. 
Actually, Ith plays a very similar role as Ic in our model. 
There is a difference however, which is that in our model 
not only the thermal but also the non-thermal leptons are 
allowed to gain (or loose) energy through Coulomb inter- 
actions with the protons. In the eqpair model the results 
are usually parametrized as a function of the non-thermal 
fraction Inth/lh where /h=/nth+/th and l^th is defined as in 
equation ((2)1. In our model such parametrization would be 
misleading since, as it will be shown below, a large fraction 
of the power provided as non-thermal electron acceleration 
(i.e. Znth ) actually ends up heating the lowest energy par- 
ticles of the distribution via the synchrotron and Coulomb 
boilers. 



2.3 Magnetic field 

While the total injection compactness I can be evaluated us- 
ing size estimates and the observed luminosity of the sources, 
the strength of the tangled magnetic field B is the most un- 
certain parameter of the problem. It is parametrized through 
the usual magnetic compactness: 



2^"o~' 



or equivalently 



(4) 



B = 5.9 X 10" 



Ib 



30Rg 20Me 



R 



M 



1/2 



(5) 



(6) 
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Figure 1. Synchrotron self-Compton models with pure non-thermal injection. Effects of varying /nth at constant inth/'s- Simulated 
photon spectra (left panel) and particle distributions for /nth= 0.3 (pink), 3 (green), 30 (red), 300 (blue), 3000 (cyan) and Inth/lB = 3. 
The other parameters are Finj = 3, 7min = li 7max = 1000, r; = 2, i? = 5 X 10^ cm. In the right hand side panel the full curves show the 
electron, while the dashed ones show the positron distributions. The dotted curve shows the shape of the injected electron distribution. 



The results of this paper will be discussed mostly in two 
different situations, namely for magnetic field in equiparti- 
tion with the radiation field and with the proton energy 
density. 

If the magnetic field is at equipartition with the radia- 
tive energy density then: 



Ur = Ltcsc/V, 



(7) 



where tesc is the photon escape time and V is the volume 
of the sphere. We can adopt the prescription of Lightman 
& Zdziarski (1987) which is used in our code and which, 
in the limit of low energy photons [hu < rrieC^), gives 
tesc — (1 + tt/3)R/c. One obtains an equipartition mag- 
netic compactness of: 



(8) 



Alternatively, if the magnetic field is in equipartition 
with the ion energy density then: 



mi 



_ 3 _ 



-to 



(9) 



where rrii/mo — 1836 is the proton to electron mass ratio, 
6i — kTi/m-iC? is the proton temperature in units of the pro- 
ton rest mass energy. fcTo = 10 MeV is a typical value for 
the temperature in the innermost part of an advection dom- 
inated accretion fiow. This gives an estimate of the equipar- 
tition magnetic field: 



Bp 



3.2 X WG 



30i?G 2OM0 Ti 



R 



M To 



1/2 



(10) 



It is also interesting to express the equipartition mag- 
netic compactness as a function of the Coulomb compact- 
ness. For magnetic field in equipartition with the proton 
energy density, we can combine equation and (|9} in or- 



der to estimate the equipartition magnetic compactness as 
a function of L: 



I Bp 



3 f Ti „ k 

oy — + 



(11) 



where y = 4tt^c is the Compton parameter. 

In the limit of large Coulomb compactness (or small 
Thomson depth), Isp — Syglc/r-^- For rx = 2 and yg— 1, 
this leads to Ibp — 0.7Zc which, in absence of other sources 
of heating or soft photons, is comparable to the values ob- 
tained assuming equipartition with radiation {Ibp = 0.4Z). 
In the limit of small coulomb compactness, however, the 
equipartition magnetic compactness saturates at 32//8 ~ 0.4 
corresponding to equipartition in a one temperature plasma 
(Bp ~ 5 X 10"" G for typical parameters). 



3 SIMULATION RESULTS 

In this section we present the incidence of the various pa- 
rameters of our model (namely /nth, ^s, ^c, Ib, ti, i?, Finj, 
7min and 7max) ou the shape of the electron and photon 
distribution. In Section |3TT] we explore models in which syn- 
chrotrons photons are the only source of soft seed photons 
(/s = 0) and the protons are cold (/c = 0) and the corona 
is powered only through non-thermal acceleration. In Sec- 
tion |3]2] we illustrate and discuss the relative importance of 
e-e Coulomb and Synchrotron self-absorption as thermalis- 
ing mechanisms. Then in Section [3. 31 we will investigate the 
effects of additional external soft-photons. Finally in Sec- 
tion [X?] we illustrate the effects of Coulomb heating by hot 
protons. The parameters and the results of the simulations 
discussed in this section are summarised in Table [1] 
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Table 1. Parameters of the model shown in Fig. [T]to[9l as indicated in the first column. The next 9 columns give the model parameters 
as described in the text. The blackbody temperature fcTbb is given in eV, R in cm. a is the resulting spectral slope in the 2-10 keV 
range. -Epeak is the maximum of the vF^ spectral energy distribution (in keV). rx is the total optical depth (i.e. including pairs). kTc is 
the temperature of the thermal component of the electron distribution (in KeV) as determined from a fit with a Maxwellian. fcT; is the 
temperature of the protons in MeV. < 7 > is the average Lorentz factor of the electron distribution. 
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Figure 2. Synchrotron self-Compton models with pure non-thermal injection. Effects of varying Ib (panel a), Finj (panel b), tj (panel 
c) and R (panel d) on the spectral energy distribution. In all panels, the result of our fiducial model obtained for (nth = 30, Ib = 10, 
Tinj = 3, 7min = Ij 7max = 1000, Ti = 2, i? = 5 X 10^ Cm is shown in red. Keeping all the other parameters constant, panel (a) shows 
the results for Ig = 1 (green) and Zs=100 (blue). Panel (b) shows the spectra for rinj=2.5 (green) and 3.5 (blue). Panel (c) presents the 
results for t; = 5 X 10~^ (pink), 0.5 (green) and 5 (blue). Finally panel (d) displays the spectra for R = 5 X 10^ (green), 5 X 10* (blue), 
and 5 X 10^^ cm (cyan). 



3.1 Synchrotron self-Compton models with pure 
non-thermal injection 

In this section we assume that the protons are cold {Ic = 0) 
and the external photons are neglected {Is = 0). Figure [T] 
shows the dependence of the photon spectrum on Znth for 
Ib = ^nth/3, which corresponds to approximate equipartition 
of magnetic field with radiation. The other fixed parameters 
are Finj = 3, 7i = 5 x 10^ cm, and ri = 2. For a wide 
range of compactness the spectrum peaks around 65 keV 
and and the X-photon index is F ~ 1.7 (see also Table [T|. 
The right panel of Fig [1] shows the steady state particle dis- 
tribution corresponding to the these photon spectra. These 
particle distributions are all thermal at low energies with 
a non-thermal tail at high energies. The non-thermal elec- 
trons produce self-absorbed synchrotron radiation peaking 
at the synchrotron turn-over frequency at photon energies 
of a few eV. The hybrid thermal non-thermal distribution 
then Compton up-scatters the self-absorbed synchrotron ra- 
diation, forming spectra extending into the hard X-ray and 
7-ray domain. These spectra are dominated by the Comp- 
tonisation by the thermalised particles. The peak and energy 
of the cut-off depends on their temperature. In the follow- 
ing we will refer to this process as Synchrotron Self-Compton 
emission (or SSC). 



In terms of particle kinetics, those results can be un- 
derstood as follows. The injected high energy electrons cool 
down by emitting SSC radiation and by interacting through 
Coulomb collisions with lower energy particles. We define 
7i as the Lorentz factor corresponding to the energy of 
the electrons radiating most of their synchrotron luminos- 
ity around the synchrotron turn over frequency. Below the 
critical Lorentz factor 7t, the particles emit synchrotron in 
the self-absorbed regime: their emission is immediately ab- 
sorbed by lower energy electrons. This very fast process of 
synchrotron emission and self absorption represents a very 
efficient and fast thermalising mechanism. Very quickly the 
low energy part of the electron distribution forms a quasi 
Maxwellian. This synchrotron boiler effect is able to transfer 
a large fraction of the non-thermal energy of the electrons 
(of Lorentz factor < 7t) into heating of the lowest energy 
particles, keeping the effective temperature around 40 keV. 

An analytical estimate of yt is derived in Appendix [K\ 
This gives 7t ~ 10, the exact value being rather insensitive 
to the model parameters, see equations (IA16|) and ([AlSp . 
Not only the electrons at or above 7t do not contribute to 
the heating of the lower energy particles but their soft syn- 
chrotron radiation also provides the dominant contribution 
to the Compton cooling of the plasma. Therefore the effi- 
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Figure 3. Effects of varying 7min in pure synchrotron self-Compton models. Simulated photon spectra (left panel) and particle distri- 
butions for 7niin= 1 (red), 1.5 (black), 2 (pink), 3 (green), 5 (blue) and 10 (cyan). The other parameters are i^th = 30, Ib = 10, Finj = 3, 
7max = 1000, T; = 2, i? = 5 X lo'^ Cm. In the right hand side panel the full curves show the electron while the dashed ones show the 
positron distributions. The dotted curve shows the shape of the injected electron distribution. 



ciency of the synchrotron boiler as a heating mechanism will 
depend strongly on the relative fraction of electrons present 
above and below 7t. At the same time Coulomb collisions 
between high energy and low energy electrons also plays an 
important (and possibly dominant) thermalising role, sim- 
ilar to that of the synchrotron boiler. The relative impor- 
tance of synchrotron and Coulomb boilers will be discussed 
in Section [33] 

At fixed Is/lnth both synchrotron and Compton losses 
scale like l^th (see equation IA20|I and the injected power 
also scales like /nth- The spectra then depend only weakly 
on compactness. At low compactness the small differences 
between the spectra of various compactness are entirely due 
to Coulomb effects and will be discussed in Section [3.21 At 
very large /nth, pair production effects become important 
and electron positron pairs contribute significantly to the 
total optical depth (tt — 9 for l^th = 3 x 10'^). As a con- 
sequence the effective temperature decreases, as more par- 
ticles share the boiler energy, but the Compton parameter 
actually increases producing even harder spectra (r=1.39 
for Znth = 3 X 10^). We also note that, at high compactness, 
pair annihilation leads to the formation of the line apparent 
around 511 keV in some spectra of Fig. [T] 

These simulations show that even though the electron 
injection is purely non-thermal, most of the radiated lumi- 
nosity emerges as quasi-thermal Comptonisation. The basic 
spectral properties of black hole binaries in the LHS appear 
to be produced over a wide range of luminosity (at least 
up to 0.3Le) through pure non-thermal electron injection in 
a magnetised plasma, with no need for additional heating 
mechanism nor soft photons. 

Let us now focus on the fiducial case Inth = 30 (L ~ 
0.03I/b) and evaluate the effects of the other parameters. 
As shown in panel (a) of Fig [2 increasing Ib increases the 
soft synchrotron emission. As a consequence, the competing 
process, i.e. inverse Compton emission of the non-thermal 



electrons, is reduced. Simultaneously, the temperature of 
the thermal component is reduced due to the enhanced syn- 
chrotron soft cooling fiux. This leads to softer X-ray spectra, 
peaking at lower energy. In the limit of vanishing Ib the elec- 
trons cannot radiate efficiently their energy through SSC. 
Other processes such as bremsstrahlung (which is otherwise 
negligible in the situations considered in this paper) and its 
Comptonisation may become the dominant cooling mecha- 
nism of the plasma. The Ib = temperature represents the 
maximum temperature achievable with pure non-thermal in- 
jection. 

Increasing Finj decreases the fraction of electrons with 
Lorentz factor larger than 7t, therefore the self-absorbed 
synchrotron emission is reduced and the energy available 
for the synchrotron boiler is increased. The temperature of 
the Maxwellian component is then higher. This leads to 
harder Comptonised spectra peaking at higher energy, with 
a weaker non-thermal fraction, as shown on panel (b) of 
Fig.[l 

Increasing n decreases the average energy of the ther- 
mal component (since more particles have to share the same 
amount of energy), this reduces the peak energy of the 
Comptonised spectrum (see panel (c) of Fig. [2}. Simultane- 
ously, at larger n the radiation energy density is larger and, 
as a consequence, the Compton losses of the non-thermal 
leptons are also larger at the expense of the synchrotron 
emission. The synchrotron self-absorbed emission is there- 
fore weaker (see equation (|A20|| ) and, in order to radiate the 
same luminosity, the Comptonised X-ray spectrum must be 
harder. 

As shown on panel (d) of Fig. [2] varying the size of 
the emitting region R has very little effects on the pho- 
ton spectrum, at least at high energy. The main effect 
is to shift the cyclotron and the self absorption frequen- 
cies. The self absorption frequency et scales like R~'' with 
b — j^Q^^p"^ see equations (|A2[I and (|A15[) . The soft cooling 
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Figure 4. Synchrotron self-Compton models with pure non-thermal injection: effects of neglecting e-e Coulomb. Simulated photon 
spectra (left panel) and lepton distributions (right panel) for Znth= 0-3 (pink), 30 (red), 3000 (cyan) and Znth/'s = 3 (same parameters 
as in FiglTJ. The full curves show the results of the full calculation. The dashed curves show the results when e-e Coulomb collisions are 
switched-off. 



synchrotron photon compactness however remains nearly 
constant (/sync oc R'^ with c = ]^o+2r"^ - ' equation (IA20|) '). 
which ensures the temperature of the thermal component 
does not change significantly. When R increases, the spec- 
trum has to be slightly softer for energy conservation rea- 
sons. Overall the effects are almost negligible except for huge 
variations of R. Panel (d) of Fig. [2] shows a spectrum com- 
puted for a size relevant to AGN {R = 5 x 10^^ cm). The 
spectrum is softer and peaks at slightly lower energy but 
overall the shape is very similar. This kind of model would 
produce LHS like spectra in AGN. 

Fig.[3]shows the effects of varying the minimum Lorentz 
factor of the injected electrons. As long as 'jmin < 2 the ex- 
act value has virtually no effects on the steady state particle 
and photon distribution. For 7min > 2, on the contrary, the 
synchrotron emission increcises quickly with '^mini cooling 
down the plasma and producing softer hard X-ray spectra. 
This effect is simply due to the fact that we keep Inth con- 
stant when increasing 7min. Since less power is injected at 
low energies, this involves increasing the normalisation of 
the injected distribution and therefore increasing the non- 
thermal particle density at energies above 7t, and therefore 
the synchrotron luminosity. This increase in normalisation is 
not linear with 7min. For "fmin = 2 the normalisation of the 
injected distribution is only about 30 percent larger than for 
7min ~ 1, and the effect is weak. For larger 7min the normal- 
isation has to increase faster. For 7min=4 and 10 it is larger 
by factors of about 2 and 5 respectively. The increased syn- 
chrotron emission has then a strong impact on the energy 
balance in the corona. 

3.2 Synchrotron boiler vs e-e Coulomb boiler 

Since the effects of synchrotron absorption on the electron 
distribution have been routinely neglected in previous works, 
it is interesting to check their importance and compare them 



with other thermalising processes. Another important ther- 
malising mechanism is the e-e Coulomb interaction which 
may keep the particle distribution close to a Maxwellian at 
low energies and transfer a fraction of the power from high 
to lower particle energies. 

The collisions between electrons/positrons and ions also 
contribute to the thermalisation process, however the effects 
are important only when the temperature of the ions is larger 
than that of the electrons (this case will be considered in Sec- 
tion [3]4]). If the temperature of ions and electrons is compa- 
rable we found the effects completely negligible compared to 
e-e collisions (and the e-p collisions were actually neglected 
in the simulations with Ic ~ presented in this paper). 

In principle bremsstrahlung emission and its self- 
absorption could also play a role similar to that of the 
synchrotron boiler. However, at least in the case of e- 
p bremsstrahlung, we found that the thermalising effects 
are fully negligible in the parameter regime considered in 
this work and we decided not take e-p bremsstrahlung 
into account in the simulations presented in this paper, e- 
e bremsstrahlung is not taken into account in our code but 
its effects are also expected to be weak (see discussion in 
Belmont et al., 2008b). 

Therefore the only important thermalising mechanism 
besides synchrotron self-absorption is e-e Coulomb. In terms 
of energy loss by non-thermal electrons the ratio of Coulomb 
losses to synchrotron losses 7c /73 oc {tt /lB)P~^'y~'^ (see Ap- 
pendix [XJ. Therefore Coulomb collisions dominate over syn- 
chrotron losses at low magnetic compactness and low parti- 
cle energies. And for this reason unless the ratio /b/tt is very 
large, the Coulomb process is the dominant process keep- 
ing the low energy part of the lepton distribution thermal. 
The synchrotron boiler is the dominant thermaliser only for 
Ib/tt > 100 (see discussion in GHS98). However since syn- 
chrotron losses dominate at higher particle energies, the syn- 
chrotron boiler can be more effective at transferring the en- 
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Figure 5. Comparison of tlie effects of the synchrotron and e-e Coulomb boilers in a SSC model with pure non-thermal injection. 
Simulated photon spectra (left panel) and electron distributions (right panel) for the full calculation (full curve) , neglecting the effects of 
synchrotron heating on the particle distribution (dot-dash), neglecting both synchrotron heating and cooling (dotted curve), neglecting 
e-e Coulomb (dashed curve) and neglecting simultaneously Coulomb, synchrotron heating and cooling (3 dots-dash). The parameters are 
those of our fiducial model (i.e. (nth = 30, Ig = 10, Finj = 3, tj = 2, = 5 X 10^ cm). 



ergy from the non-thermal to the thermal particles (i.e. as 
a heating mechanism). In practice, the synchrotron boiler 
contributes significantly to the heating of the Maxwellian 
distribution for Ib/tt > 0.2 (see equation HA22|) ). 

Using the results from our code, we can estimate the 
synchrotron and Coulomb heating power. We define it as 
the net power transferred from particles of Lorentz factor 
above 7 = l-f lOkTc/rricC^ to particles of lower energy. Con- 
sidering the models of Fig. [T] we find that for Ib =0.1 the 
synchrotron heating is negligible compared to Coulomb (in 
terms of compactness the Coulomb and synchrotron heating 
are 0.19, and 3xl0~^ respectively). The model with Z_b = 1 
is just above the threshold for which synchrotron heating 
becomes significant. For this model, we find 1.48 and 0.14 
respectively for the Coulomb and synchrotron heating com- 
pactness. In the model with Ib = W they are comparable 
(5.9 and 4 for Coulomb and synchrotron heating respec- 
tively). At /_B = 100, heating is dominated by synchrotron 
(we find 26 and 59 for the Coulomb and synchrotron heating 
compactness respectively). Finally for Ib ~ 1000, Coulomb 
effects are very weak (we find 44 and 314 for the Coulomb 
and synchrotron heating compactness respectively). 

Besides thermal heating, Coulomb interactions also 
have a significant effect on the non-thermal tail of the lep- 
ton energy distribution. Indeed, for the non-thermal parti- 
cles of lowest energy. Coulomb cooling usually dominates 
over Compton and synchrotron losses. This tends to deplete 
faster the electrons from this energy range. As a result the 
electron distribution is less steep than the expected Finj + l 
slope obtained in the case of pure radiative losses (see equa- 
tion (|A12|) ). This flattening occurring in the range p =1-10, 
can be observed in the electron distributions presented in 
Fig. [1] This effect is also illustrated by Fig. |4] which com- 
pares particles and photons spectra of the simulations with 



/s=0.1, 10 and 1000 of Fig. [T]with models computed with 
the e-e Coulomb collisions switched-off. 

Because the importance of Coulomb losses decreases 
with compactness, the deviation from a power law starts at 
higher electron energies for lower compactness. An estimate 
for the momentum corresponding to the transition from ra- 
diative to Coulomb cooling is provided by equation HA9[) . For 
the simulations with Z_b=0.1, 10 and 1000 this gives pc=19.4, 
1.98 and 0.56 respectively. 

The main effects of Coulomb collisions on the photon 
spectrum is to increase the sharpness of the thermal comp- 
tonisation cut-off by depleting the supra-thermal particles. 
This effect is very significant, and should be taken into ac- 
count when comparing to data. For instance, in section W7I\ 
we compare SSC models with a LHS spectrum of Cygnus 
X-1 and find a good agreement when Coulomb cooling is in- 
cluded while no reasonable fit could be found neglecting this 
effect. As seen on Fig. 2] the effects of Coulomb collisions 
on the spectral shape decrease at higher compactness as dis- 
cussed above. For the high compactness case Ib = 1000, the 
effects of Coulomb collisions are practically negligible, both 
thermalisation and thermal heating are dominated by the 
synchrotron boiler. 

In addition, as expected, the temperature of the thermal 
electrons increases when Coulomb is taken into account (see 
Table [T|. As a consequence, the vFi, spectra peak at slightly 
higher energies in the simulations including Coulomb. The 
decrease of the temperature with compactness in the simu- 
lations of Fig [T] is related to the dependence of Coulomb 
heating on Ib- In the simulations without Coulomb, the 
lower temperature of the Maxwellian is compensated by the 
much larger density of supra-thermal particles which also 
contributes to the Comptonisation process. As the average 
energy of the Comptonising particles is almost identical. 
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Figure 6. Non-thermal injection models dominated by external photons. Effects of varying the non-thermal compactness l^tb a^t constant 
Is = 15 and Ib = 10. on the photon spectra (left hand side) and the lepton distributions (right hand side). The pink, green and blue 
curves stand for Znth=5, 15 and 50 respectively. The full curves are models with Thomson depth Ti=2, while the dot-dashed curves show 
the results for Ti=0.2. The other parameters are fcT^ij = 420 eV, Finj = 3, 7niin = Ij 7max = 1000, R = 5 X 10^ cm, see Table [Tl 



the X-ray spectral slopes are not strongly affected by the 
Coulomb collisions. 

Moreover, we note that the Coulomb and synchrotron 
boilers are competing to produce the same result, when 
Coulomb collisions are neglected the synchrotron boiler 
takes over transferring a large fraction of the non-thermal 
energy that would have otherwise been transferred via 
Coulomb collisions to the thermal pool of electrons. In the 
case, Ib = 0.1, the synchrotron heating power is 7.2 x 10~^, 
i.e. larger by a factor 24, when the Coulomb boiler is ne- 
glected. In our fiducial model Ib ~ W the synchrotron heat- 
ing power doubles to become 7.4, compensating almost fully 
for the lack of Coulomb heating (the total heating compact- 
ness was about 10 in the full calculation). 

In order to understand better the respective roles of 
Coulomb and synchrotron in the thermalisation and plasma 
heating processes, it is useful to see what happens when the 
synchrotron boiler is turned off. FigOshows our fiducial sim- 
ulation for Znth=30 in which the heating of the leptons due to 
synchrotron self-absorption is neglected (but heating by e-e 
collisions is taken into account). In this simulation, the en- 
ergy of the absorbed synchrotron photons is artificially lost 
and a large fraction of the injected energy is not radiated. 
As a consequence, the luminosity is lower by a factor of 4 
with respect to the full calculation. The temperature of the 
thermal leptons is then significantly lower (only 16 keV see 
Table[T}. As a consequence the spectra is much steeper, with 
photon index F ~ 2, incompatible with the typical slopes ob- 
served in the LHS. The heating of the thermal population 
by e — e Coulomb collisions appears too weak to sustain the 
electron temperature observed in the hard state. 

However, when comparing to Coulomb process we 
should consider the net effect of the magnetic field on the 
lepton distribution and switch-off both synchrotron heat- 
ing and cooling. Fig [5] shows the same simulation with not 
only the synchrotron heating but also synchrotron cooling 



turned-off. Then the temperature is higher (fcTe=26 keV) 
than in the simulation including synchrotron cooling yet sig- 
nificantly lower than in the full calculation. However this is 
essentially because the synchrotron losses of the high energy 
electrons are not taken into account and therefore the den- 
sity of electrons above 7t is larger by about a factor of 2. 
As a consequence, the self-absorbed synchrotron luminos- 
ity is also larger by a factor of two. In this simulation the 
thermal Coulomb heating compactness is 8.6 i.e. compara- 
ble to total (Coulomb -f- synchrotron) heating rate of the 
full calculation (~ 10). What makes the temperature lower 
is the larger Compton cooling rate of the thermal electrons. 
This leads to a softer hard X-ray photon index of F ~ 1.8. In 
fact the temperature of the thermal electrons is close to that 
obtained when neglecting e-e Coulomb, however, because of 
the enhanced Compton cooling rate the average Lorentz fac- 
tor is lower ({7) = 1.11 when synchrotron heating and cooling 
are neglected, {7) = 1.13 when Coulomb is neglected). 

Finally, we note that the results without Coulomb are 
close to those of the full calculation. On the contrary, the 
photon spectra obtained neglecting both synchrotron cool- 
ing and heating are similar to those obtained neglecting all 
thermalising processes (i.e. switching-off both Coulomb and 
synchrotron boilers also shown in Fig. [5]). 

3.3 Effects of external soft photons 

In addition to soft photons produced internally through syn- 
chrotron, the hot comptonising plasma may intercept a frac- 
tion of the thermal radiation from the accretion disc. We 
denote as Is the compactness associated to this additional 
soft photon injection. 

The main effect of the external soft photons field, is 
to increase the Compton cooling rate of the leptons. The 
equilibrium temperature of the thermal component is lower 
than in pure SSC models. If k is strong (Zg ~ ^nth) the ther- 
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Figure 7. Non-thermal injection models dominated by external soft photons. Effects of varying the magnetic compactness Ib at constant 
is=15 and Znth=15i on the photon spectra (left hand side) and on the lepton distributions (right hand side). The pink, green and blue 
curves stand for is=0, 10 and 50 respectively. The full curves are models with Thomson depth Ti=2, while the dot-dashed curves show 
the results for Ti=0.2. The other parameters are fcT^ij = 420 eV, Finj = 3, 7min = Ij 7max = 1000, R = 5 X 10^ cm, see Table [Tl 



malised leptons are so cool that most of the luminosity is 
radiated by the non-thermal particles. Basically the effect 
of Is, is to increase the fraction of radiation emitted through 
non-thermal Comptonisation. 

At fixed la and Ib increasing /nth increases the amount 
of energy available for the synchrotron boiler, the tempera- 
ture of the plasma increases and the high energy emission is 
not only stronger, but also harder as can be seen in Fig. |6l 
As shown on this figure and in Table [TJ the temperature is 
higher at lower optical depth because less particles share the 
same energy. This behaviour is similar to what obtained in 
the SSC model. In the SSC model this resulted in a softer 
spectrum at low optical depth (see Fig. ^) while Fig. [S] 
shows harder X-ray spectra at lower optical depth. This is 
because, as discussed in section I3TT1 in the SSC models the 
flux of seed photons for comptonisation is larger at small rx, 
whereas it is independent of tt when Is is important. 

We note that in presence of soft photons the effects 
of the magnetic field are somewhat different from those 
of pure synchrotron self-Compton models. As discussed in 
Section 13.11 and illustrated in Fig. [2^, in absence of soft 
photons, increasing Ib increases the radiative losses (both 
trough synchrotron and and synchrotron self-Compton), and 
causes the temperature of the plasma to decreases. When the 
losses are dominated by Compton on external soft-photons, 
then increasing Ib has negligible effects on the radiative 
losses of the plasma, yet this transfers a larger fraction of 
non-thermal energy to the lower energy particles which re- 
sults in a higher plasma temperature. Therefore, when syn- 
chrotron self-Compton dominates, the magnetic fields tends 
to cool down the thermal component, when external Comp- 
ton dominates it heats up this distribution. This effect how- 
ever strongly depends on the Thomson optical depth of the 
plasma. For instance if we consider the simulation shown 
in Fig. [7]with Ib = 10, Za=15, and /nth=15 the equilibrium 
temperature of the thermal component is about 14.6 keV. 



Setting Zs = 0, we obtain the same temperature, while set- 
ting Ib = 50 we obtain kTe = 15.3 keV, i.e. the temperature 
increases by only 5 percent. Then, repeating the same ex- 
periment with an optical depth reduced to 0.2, we find a 
temperature of kTa—SG keV for Zs = 10 and the temper- 
ature increases by a factor of 2 between the case Z_b = 
(fcTc=59 keV) and Zi3=50 (fcTe^ll? keV). This is because, 
while the Compton and synchrotron losses are nearly inde- 
pendent of the optical depth, the energy of the synchrotron 
boiler is shared among a smaller number of particles in the 
optically thin case, making the synchrotron heating much 
more effective. We also note that even when the optical 
depth is too thick for the magnetic field to play a siginfi- 
cant role in heating the thermal component of the plasma, 
it significantly affects the emission of the non-thermal parti- 
cles. For a stronger Ib the high energy particles radiate their 
energy in the UV band rather than at gamma-ray energies 
(see Fig. [7]). Therefore increasing Ib also reduces the frac- 
tion of high energy radiation emitted through non-thermal 
Comptonisation . 

Fig [8] shows a sequence of spectra, starting from our 
fiducial case Znth ~ 30, Ib = 10 in which Is is increased 
keeping the total compactness Z = Zs -I- Znth constant. As 
the soft photons flux rises in the corona, inverse Compton 
increasingly dominates the cooling of the non-thermal par- 
ticles eventually turning the synchrotron boiler off. As Is is 
increased, the electron distribution cools down to the Comp- 
ton temperature. The emissivity of the thermal component 
is strongly reduced and the hard X-ray spectrum becomes 
gradually dominated by the emission of the cooling non- 
thermal distributiorlfl. Then, the non-thermal emission van- 
ishes as Znth decreases. We note that this result as well as the 



^ The secondary bumps observed in the electron distribution of 
some of the simulations presented in Fig|6l [T] |8] are essentially 
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Figure 8. Pure non-thermal models with external soft photons. Effects of varying the soft photon compactness Is at constant I on the 
photon spectrum (left hand side) and the particle distributions (right hand side). The red curve shows our fiducial model obtained for 
= 0, ^nth = 30, Ib = 10, Finj = 3, 7min = li 7inax = 1000, T; = 2, i? = 5 X lo'' Cm. The other curves show a sequence of models obtained 
keeping I = Ib + Inth = 30 while the soft photon compactness Ib is varied. The black, pink, green, blue, cyan and purple curves show the 
results for Ib=3, 9, 15 21, 27, 29.7 respectively. The soft photons are injected with a blackbody spectrum of temperature fcTijb=213.5 
ll^* eV. In the right hand side panel, the full curves show the electron while the dashed ones show the positron distributions. 



main results of Section [3. II are qualitatively similar to those 
obtained independently by Vurm and Poutanen (2008). 

This spectral sequence, where Is is increased, is remi- 
niscent of that observed during the state transition of ac- 
creting black holes (DGK07; Del Santo et al. 2008). This 
indicates that cooling by the soft photons from the accre- 
tion disc probably plays a major role in the changing ap- 
pearance of the corona during such events. It was known 
from spectral fits with eqpair that due to the strong ther- 
mal emission from the disc, the ratio Ih/L was lower in the 
HSS, implying a lower temperature of the plasma. Here we 
show that the higher Inth/lh ratio of the HSS is naturally 
explained with this synchrotron boiler model. However as 
will be discussed in Section [5] real spectral state transitions 
in X-ray binaries are a bit more complex and the spectral 
changes cannot be entirely attributed to a change in the soft 
photon compactness. 



3.4 Effects of e-p Coulomb heating 

In this section we investigate the effects of Coulomb heating 
by a thermal distribution of protons in a two-temperature 
plasma. Fig [9] shows the evolution of the spectrum when Ic 
is varied at fixed I = l^th + 4 = 30, Ib = 1/3, la = 0. The 
other parameters are those of our fiducial model. When the 
fraction of Coulomb heating increases, the effects of the non- 
thermal component both in terms of high energy emission 
and soft cooling photon flux are reduced. As a consequence 
the temperature of the thermal component increases. This 



due to the use of the p'r(p) representation in which a power law 
distributions in 7 space has a maximum around p = 1. 



leads to an increasingly hard spectrum peaking at increas- 
ingly large photon energies. Finally, when Coulomb heating 
dominates, the equilibrium electron distribution is very close 
to a Maxwellian and, as a result, the non-thermal high en- 
ergy tail is absent from the Comptonisation spectrum. Due 
to the effects of the pair production thermostat (Svensson 
1984), the temperature of the electron plasma is limited to 
a few hundred keV. The temperature of the hot protons 
(shown in Table[TJ increases with Ic in accordance with equa- 
tion ((4|. 



4 APPLICATION TO CYGNUS X-1 

In this section, we compare the results of our simulations 
to high energy spectra of Cygnus X-1. Since its discovery 
in 1964 (Bowyer et al. 1965), this bright persistent source 
has been intensively observed by all the high energy instru- 
ments, from soft X-rays to 7-rays. It was the first dynami- 
cally proven black hole, and still remains the prototype of 
black hole candidates. The mass of the black hole is in the 
range 14-29 M© (Ziolkowski 2005) and its distance is 2.15± 
0.2 kpc (Massey, Johnson, & Degioia-Eastwood 1995, see 
also the discussion and references in Ziolkowski 2005). 

In most of the following we assume a black hole mass 
of 20 A'Iq at a distance of 2 kpc. Most of the X-ray emission 
is expected to be produced in a region of size R ~ 30i?G — 
10* cm, where most of the accretion power is released. This 
is what we will assume in our simulations. However the ob- 
served luminosity may also arise from several distinct and 
more compact active regions (as in the patchy corona model 
for instance) . In order to take into account the uncertainties 
in the source parameters and geometry, we denote: 
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Figure 9. SSC models with heating by hot protons. Effects of varying the Coulomb heating compactness /c at constant / on the photon 
spectra (left hand side) and the particle distributions (right hand side). The red curves show our fiducial model obtained for Ic = 0, 
'nth = 30, Ib = 10, Finj = 3, 7inin = Ij 7max = 1000, T; = 2, R = 5 X 10^ Cm. The other curves show a sequence of models obtained 
keeping I = Ic + 'nth = 30 while the Coulomb heating compactness Ic is varied. The black, pink, green, blue and cyan curves show the 
results for 1^=10, 15, 20, 29, 30 respectively. In the right hand side panel, the full curves show the electron distributions while the dashed 
ones show the positron distributions. 



5 = 



30-Rg 20Mq d? 
nR M dl ' 



(12) 



where n is the number of distinct active regions in the ac- 
cretion flow [R is their typical size), and d is the distance to 
the source and do=2 kpc. Then for a given observed flux F, 
the compactness parameter can be estimated as; 



I ■ 



O 



(13) 



where Fo=6.8xl0~^ erg s~^ cm~^. 

In order to test the synchrotron boiler model against 
spectra of accreting black holes it is important to have good 
quality data extending above a few hundred keV and up 
to at least several MeV in order to estimate the flux level 
of the non-thermal tail. With this regard, the average hard 
state spectrum of Cygnus X-1 obtained with the Compton 
Gamma-ray Observatory (CGRO) probably represents the 
best data so far for a hard state source at MeV energies. We 
use this CGRO spectrum combined with a typical Beppo- 
SAX hard state spectrum in order to extend the energy cov- 
erage down to the soft X-rays. Similarly we will also compare 
our results to the combined BeppoSAX / CGRO spectrum of 
the 1996 soft state of Cygnus X-1. These two spectra are 
shown on Fig. 1101 Both sets of data were extensively pre- 
sented and analysed in MC02. These authors fit both spec- 
tra with the hybrid thermal-non-thermal models eqpair and 
found a very good agreement. 

However, as discussed earlier, the eqpair model does 
not take into account the effects of magnetic field nor the 
Coulomb interaction with protons. In order to constrain the 
magnetic field and proton temperature we would like to fit 
those spectra with the results of our code. However our code 
has not been optimised yet and is too time-consuming in or- 
der to perform detailed spectral fitting. As a preliminary at- 



tempt to compare SSC models to spectra of accreting black 
holes, we will adopt a simpler approach and compare only 
qualitatively the shape of the simulated spectra with the un- 
folded data. Obviously, it will not be possible to compare the 
'goodness' of different models nor evaluate the uncertainties 
on the model parameters. Nevertheless we will see that these 
simple 'fits by eyes' already provide interesting constraints. 
For both spectra, we will proceed in three steps: 

(i) First, we use the best fit eqpair parameters found by 
MC02, replacing the unspecified thermal heating of eqpair 
by Coulomb heating (i.e. setting Ic = /th and Ib = 0), we 
check that our code is able to reproduce the shape of the 
spectrum, with a model without magnetic field. This allows 
us to estimate the proton temperature implied by the best 
fit models of MC02. This proton temperature depends on Ic 
(as can be approximated by equation Q). 

(ii) However, for a fixed (observed) X-ray flux, Ic depends 
on the size and number of emitting regions as well as on 
the distance of the source (I oc 5). In their spectral fitting 
procedure MC02 used an additional normalisation constant 
(directly related to 5) which is a free parameter account- 
ing for uncertainties on the size and distance and actually 
disconnecting the compactness parameter from observed lu- 
minosity. The best fit compactness parameters of MC02 lead 
to 5 significantly different from unity. We therefore compute 
a similar model with all compactness parameters divided by 
5 (so that the observed luminosity is produced for 5=1) 
and then check that we still have a good agreement with 
the data. If necessary, we modify the parameters to improve 
the fit and derive what we believe is a more realistic proton 
temperature. 

(iii) We then set the Coulomb heating Ic = Q and, allow- 
ing for the presence of the magnetic field, we explore the 
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Table 2. Models parameters of the simulations shown in Fig. 1111 In addition to the models parameters the resulting proton temperature 
as well as the temperature of the thermal component of the electron distribution obtained from a fit with a Maxwellian are given. The 
last two columns give the ratio of magnetic to proton energy density (computed according to equation and the magnetic to radiation 
energy density (computed using the radiation energy density resulting from the simulation). 



Ref. 


'nth 


'c 


's 


Ib 


Ti 




7min 


5-1 




fcTc (keV) 


kTi (MeV) 




Ib/Ib„ 


HI 


2.51 


28.1 


1.8 





1.34 


2.00 


1.5 


0.13 


1.45 


97 


13 








H2 


0.39 


4.36 


0.28 





1.45 


2.00 


1.5 


1 


1.45 


92 


1.6 








H3 


1.00 


3.00 


0.29 





1.45 


3.00 


1.5 


1 


1.45 


85 


1.0 








H4 





4.75 





200 


1.45 






1 


1.45 


83 


1.3 


37 


131 


H5 


4.75 








0.475 


1.45 


3.50 


1 


1 


1.45 


76 




1.5 


0.31 


SI 


0.36 


0.17 


3.20 





0.11 


2.64 


1.5 


5.5 


0.11 


76 


9.4 








S2 


1.98 


0.96 


17.6 





0.11 


2.64 


1.5 


1 


0.12 


67 


50 








S3 


5.50 





16.7 


17.6 


0.11 


2.10 


1 


1 


0.12 


48 




1129 


3.2 



possibility of reproducing the spectrum with a pure non- 
thermal acceleration model. 

For the purpose of our comparison with the data and 
in order to make our results comparable with the spectral 
fits of MC02, we now model the soft thermal emission of 
the cold accretion disc using a multicolour blackbody disc 
model (diskbb in XSPEC, Mitsuda et al. 1984) rather than 
a simple blackbody spectrum. When comparing our simu- 
lations with the observed spectra we account for the pres- 
ence of a reflection component which is calculated using the 
PEXRiv routine in XSPEC (Magdziarz & Zdziarski 1995). 
We also take into account the standard column density, 
A'h = 6 X lO'^^ cm, of absorbing material in the direction 
of the source (we use the wabs model in XSPEC). The re- 
sulting model parameters are given in Table [2] and the un- 
absorbed model spectra are displayed in Fig llll 

4.1 Low-hard state 

Cygnus X-1 spends most of the time in the LHS, with a 
rather stable flux of about 3.5x 10~* erg s"^ cm~^. This sets 
a constraint on the compactness I ~ 5S (see equation (|13|l V 
Spectral fits with thermal Comptonisation models in- 
dicate that, in the LHS, the optical depth is in the range 
tt ~ 1.5 — 3 and the electron temperature 50 to 100 keV 
(DGK07). Then equation provides an upper limit on the 
proton temperature of about 0.1+1.9(5 MeV. This temper- 
ature would be achieved if the electron heating was fully 
due to Coulomb interaction. If any other form of heating 
or acceleration is effective in Cygnus X-1 then the proton 
temperature must be lower. However, the observations show 
a non-thermal tail at high energy in excess of pure thermal 
Comptonisation model. While we have no direct evidence of 
Coulomb heating, the high energy excess indicates that non- 
thermal acceleration mechanisms are at work in the source. 
The temperature of the protons in the X-ray emitting re- 
gion of Cygnus X-1 is therefore very likely to be much lower 
than typical ADAF temperatures (10-100 MeV). Our nu- 
merical simulations suggest even lower proton temperatures 
(see Table 0. 

At equipartition with radiation, equation ((8)1, gives 
Ib^i ~ 2(5 (i.e Br ~ 8 x W^5^/n do/d G) , and if equipar- 
tition is achieved with the protons then 0.37 < Ib^ < 
0.44-1-7.9(5 depending on the importance of Coulomb heating 
(3.6 X 10^5-^^'^ < BpS-^n-^^^ d/da < 1.7 x 10® G). 



4-1.1 Best fit EQPAIR model 

For the hard state data, MC02 froze the disc temperature 
at 200 eV, and fixed 7min = 1-5, 7max=10"^. Their best fit 
parameters are then Ih/ls ~ 17^3, Znth/^h ~ 0.082^0 032 5 
Finj = 2.0to'.2' ''"i ~ l-34lo'5 s-iid the soft compactness 
Ib = 1.8l^ g. They found a refiection component with a rel- 
ative amplitude n/2Tv = 0.52lE;;g^ was required (reflection 
on neutral material was assumed). 

4-1.2 Models with e-p Coulomb heating 

We used our code to simulate a similar model (referred as 
model HI in Table [21 . We used the same parameters as 
those found by MC02 except that we replaced the unspec- 
ified 'thermal' heating compactness of eqpair by Coulomb 
heating with Ic — hh- The magnetic compactness was set to 
Ib = and therefore the thermalisation occurs through the 
e-e Coulomb boiler. We found a spectrum which is very close 
to the EQPAIR one, and that provides a good representation 
of the data. This spectrum is shown in Fig. 1111 As shown 
in Table [2] for this simulation (HI) the equilibrium proton 
temperature is 12.9 MeV. 

However, the compactness parameters of MC02 (Z —20), 
when compared to the observed fiux, imply ~ 0.13, 
which we believe is unlikely. Our preferred size and distance 
of the source (i.e. S = 1) imposes a much lower value of the 
compactness (/ ~ 5) in order to reproduce the observed flux. 
If we adopt this value (model H2 in Table [2| we obtain a 
proton temperature of only 1.6 MeV. However, in this case 
the agreement with the CGRO spectrum is not as good. The 
shape of the high energy excess is poorly reproduced. 

This is because the effect of e-e Coulomb is more im- 
portant in the lower compactness model and flattens the 
slope of the non-thermal excess (as discussed in Section r3.2p . 
Therefore, for a lower compactness an injection slope Finj 
steeper than the one of MC02, is needed to fit the data. 
For instance, a better agreement could be found by setting 
Finj = 3, Znth = 1, Zc = 3, which leads to a proton tempera- 
ture of 1 MeV (model H3 in Table [2]) . The photon spectra 
corresponding to these three models can be compared in 
Fig. [Ill 

For comparison. Fig. [TT] also shows a pure thermal SSC 
model {Is = 0, Znth = 0, model H4 in Table [5} leading to 
a proton temperature of about 1.26 MeV. The agreement 
with the data requires an excessively super equipartition 
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magnetic field (/b//sp~30) and of course does not repro- 
duces the non-thermal high energy excess. This confirm the 
resuhs of Wardzihski & Zdziarski (2000 here after WZOO) 
who derived analytical estimates of the thermal SSC lumi- 
nosity, compared them to the luminosity of several black hole 
sources and showed that the magnetic field would have to be 
implausibly strong in order to match the observed luminos- 
ity. For Cyg X-1 however, they found that for equipartition 
magnetic field and typical parameters, thermal SSC repre- 
sents a significant fraction, and could possibly dominate the 
observed luminosity. This is not in contradiction with our 
results. In fact, these authors estimate the same magnetic 
field as suggested by our results, but they assume a pro- 
ton temperature ten times larger (10 MeV) and an emitting 
region 3 times smaller, so that they find a proton energy 
density 30 times larger. Although in this case, the magnetic 
field is indeed in equipartition with matter, we have shown 
that the proton temperature must be smaller, which implies 
a magnetic field clearly above equipartition. 



4-1.3 Model with pure non-thermal acceleration 

Alternatively, the hard state spectrum can be quite well rep- 
resented with a pure SSC model which does not include any 
form of thermal or Coulomb heating. Fig 1101 shows a com- 
parison of the observed spectrum with a simulated spectrum 
with Ib = 0.475, /nth = 4.75, h = 0, h = 0, Finj = 3.5 (model 
H5 in Table [2} . A steep power-law electron injection slope 
is required in order to reproduce the hard slope of the hard 
X-ray spectrum (for flatter injection, the soft cooling syn- 
chrotron flux is too strong). It also helps fltting the steep 
slope of the non-thermal excess. The magnetic field is ap- 
proximately at equipartition with the particles and about 30 
percent of equipartition with the radiation field. Note that 
in this model we set the minimum injection Lorentz factor 
of the electron to 1 (while MC02 use 1.5). Again this has 
very weak efi'ects on the resulting spectrum. 

We also note that in this model the radiation energy 
density dominates over that of the gas, which might be a 
cause of instability. This can be avoided if the size of the 
emitting region is larger than what we assumed by at least 
a factor of 5 (e.g. 150 -Rg instead of 30 -Rg). Finally, it 
is apparent in Fig. [10] that our model underestimates the 
flux below 1 keV. This may be due to the fact that our 
model does not include any disc thermal component. Such 
a thermal component in the LHS was firmly established by 
previous studies. However, as discussed in Section [5.11 it is 
possible that this thermal disc photon flux does not intercept 
the X-ray emitting region and has no effect on the cooling 
of the coronal electrons. 



4.2 High soft state 

In the HSS, the source is about 4 times more luminous than 
in the hard state (see MC02). This leads to I ~ 205. At 
equipartition with radiation the magnetic field compactness 
is Ib^ =i 5(5 (or Br ~ 1.5 x lO'^S^/n d^/d G). 
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Figure 10. A comparison of the average CGRO spectra for the 
HSS (blue) and LHS (red) of Cygnus X-1 from MC02, with mod- 
els involving only injection of non-thermal particles as heating 
mechanism. At low energy the CGRO data are complemented by 
BeppoSAX data (see MC02 for details). The thick purple line 
shows our absorbed soft state model S3 which parameters are 
given in Table [2] while the orange line shows the hard state model 
H5 of the same table. Reflection components were added to both 
spectra and are shown by the thin dot-dashed and dashed curves 
for the soft state model and the hard state model respectively. 



4-2.1 Best fit EQPAIR model 

The soft state spectrum was fit by MC02 with eqpair with a 
disc temperature of 370 eV, and the best fit parameters were 



O.III0021 ^^'^ ^ ionised reflector (with ionisation parame- 
ter ^ = lOOlgo" a.nd reflection amplitude n/27r = 1.3lo;3). 
Because , in their fit, compactness is decoupled from the ob- 
served luminosity (i.e. 5 is a free parameter), they find the 
soft photon compactness is poorly constrained. They find 
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4-2-2 Models with Coulomb heating 

We simulated a similar model replacing Ith by Ic (model SI 
in Table [2} and again we found a spectrum that was very 
similar to that of eqpair and consistent with the Beppo- 
SAX / CGRO data. For this set of parameters we derive a 
proton temperature of 9.35 MeV. The compactness param- 
eter then implies 5~^=5-5- We note that according to the 
fit of MC02, in the soft state is larger than in the hard 
state by a factor of 42. This would imply that nR is larger 
in the soft state, which appears in conflict with the timing 
data suggesting the size of the emitting region is smaller in 
the HSS (see e.g. DGK07). For our assumed size of the emit- 
ting region of 10** cm and a distance of 2 kpc (i.e. 5 ~ 1) , 
we need a soft photon compactness Is — 17.6 to match the 
observed flux. The resulting proton temperature is then 50 
MeV (see model S2 in Table [2J. 



4-2-3 Model with 



non-thermal acceleration 



Alternatively, a qualitative agreement with the average soft 
state spectrum of MC02 can be obtained without any 
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Coulomb heating (Zc ~ 0). Fig llOl shows a comparison of the 
data with a simulated spectrum with Ib ~ 17.6, h ~ 16.7, 
^nth = 5.5, Finj — 2.1 (model S3 in Table [2]). This suggests a 
relatively strong magnetic field of 2.3 x 10® G (about 3 times 
equipartition with radiation). Lower magnetic fields do not 
give a satisfactory description of the observed spectrum be- 
cause the thermal component of the electron distribution is 
then too cool (the shape of the spectrum in the 2-20 keV 
band is not reproduced). A higher magnetic field seems re- 
quired in order to heat up the plasma (in the way discussed 
in Section r3.3|l . 

Our results for the soft state show that the spectrum 
can be reproduced either assuming cold protons and a mag- 
netic field above equipartition, or no magnetic field and a 
relatively high temperature of the protons. The magnetic 
and the proton thermal energy density of the two models 
are comparable. Obviously intermediate models with both 
magnetic field and hot protons are not excluded and would 
allow for a lower proton temperature and/or magnetic field. 



5 DISCUSSION 

The main result of this paper is that both the LHS and and 
the HSS spectra can be produced with pure non-thermal 
acceleration models with no need for thermal heating (as 
assumed in the eqpair model) or hot protons (as assumed 
in ADAF-like models). In the framework of accretion disc 
corona model, this explains how the electron distribution 
could look so different in both spectral states while the en- 
ergy dissipation processes remain the same (i.e. non-thermal 
acceleration). In the context of the hot accretion flow mod- 
els for the hard state, our result explains why the electron 
distribution is essentially thermal even if, as expected, non- 
thermal acceleration/heating processes are important (as 
also shown by Mahadevan & Quataert 1997). In the fol- 
lowing we summarise and discuss further the constraints on 
the main parameters of the model. 



5.1 Soft photon compactness 

5.1.1 Low hard state 

In the LHS of several sources, it is likely that Compton scat- 
tering off soft disc photons contributes to the high energy 
emission and could even dominates over SSC. This is sug- 
gested by the correlation that is observed in Cygnus X-1 and 
other sources between the strength of the disc reflection fea- 
tures and the spectral slope of the X-ray emission, indicating 
a stronger Compton cooling when the disc is 'closer' to the 
X-ray emitting region (Zdziarski, Lubiriski & Smith 1999; 
Gilfanov, Churazov & Revnivtsev 1999). Such disc thermal 
emission is also clearly detected in the LHS. However it is 
possible that it does not intercept the Comptonising cloud. 
This would be the case if this emission arises from a disc 
truncated at large distances from the central region where 
most of the Comptonised emission is produced. In any case 
our results demonstrate that such external soft seed photons 
are not required by the data. The broadband high energy 
spectrum of Cygnus X-1 appears to be consistent with a 
pure SSC model with a slightly sub-equipartition magnetic 



fleld. Our results, of course, do not exclude that Comptoni- 
sation of external soft photons could dominate in the LHS. 
This would simply require either a lower magnetic field, or 
some heating by hot protons in order to keep the electron 
temperature around 100 keV. 

5.1.2 High soft state 

We have shown that in the HSS the presence of a strong 
illumination by soft photons has important consequences on 
the coronal electron distribution. This soft photon flux cools 
down the thermal particles so that most of the hard X-ray 
luminosity is produced by the non-thermal electrons. 

5.1.3 Spectral transition 

The main parameter controlling the spectral evolution dur- 
ing state transitions is the soft photon compactness. The 
change from quasi-thermal to essentially non-thermal Comp- 
tonisation spectrum is achieved naturally once the effects of 
synchrotron boiler are taken into account. However, the data 
indicate that other coronal parameters change during spec- 
tral state transitions. Broad band spectral fits of HSS spec- 
tra of Cygnus X-1 with eqpair by G99 and MC02 suggest 
the Thomson optical depths is lower by at least a factor of 
5 in the HSS. The situation is therefore more complex. Our 
comparison with the average soft and hard state spectra of 
Cygnus X-1 (presented in Section |3| also shows that other 
coronal parameters must be changed in order to reproduce 
the details of the observed spectra. 

5.2 Magnetic field 

It is interesting to express the magnetic field values we found 
in our comparison with the data as a fraction of the equipar- 
tition compactness (radiation or particles). Indeed, because 
the shape of the spectrum depends on the ratio Ib/1 rather 
than on the absolute value Ib, the ratios and Ib/Ib^i is inde- 
pendent of the uncertainties on the source size and distance 
(i.e. 8) and Ib/Ibp depends on 5 only at low proton temper- 
ature. 

5.2.1 Low hard state 

We qualitatively reproduce the LHS spectrum of Cygnus X-1 
with a pure non-thermal model in which the magnetic field is 
close to equipartition with the particle energy {Ib — 3y/8). 
We note that Wardzihski and Zdziarski (2001) (hereafter 
WZOl) using analytical estimates for the SSC luminosity 
accounting of the synchrotron emission of the non-thermal 
electrons, found that the magnetic field in Cyg X-1 must 
be strongly sub-equipartition with the protonQ. This is be- 
cause they assumed a much higher proton temperature of 
10 MeV, while in our model without Coulomb heating (H3) 
the protons have the same temperature as the electrons (80 

* On the contrary WZOO assumed the electron distribution is a 
pure Maxwellian. They found that the magnetic field would then 
be large, which is in agreement with our result for a pure Coulomb 
heating model, as discussed in section r4.1.2l 
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Figure 11. Various models providing a qualitative agreement with the observed high energy spectrum of Cygnus X-1, as described in 
Table [2] The left hand-side panel shows the spectra while the right hand side panel shows the lepton distributions. The models for the 
soft state are shown by the thin curves (model SI: dot-dashed (red); S2: dashed (green); S3: full (purple)). The models for the hard state 
are shown by the thick curves (model HI: thick dashed (cyan); H2: dot-dashed (pink); H3:dotted (blue); H4: triple-dot-dashed (black); 
H5: full (orange)) 



keV). More generally the observed luminosity sets an up- 
per limit to the proton temperature of about 25 MeV as 
shown above. Actually, our constraints are fully consistent 
with those of WZOl. In any case, the magnetic field infered 
from model H3 is probably an upper limit on the actual 
magnetic field in the source. The presence of external soft 
photons would imply a lower B as mentioned just above. If 
Coulomb heating is important then /nth must be lower to 
conserve the same luminosity. At the same time, in order 
to keep the same non-thermal Compton luminosity it would 
be necessary to reduce the synchrotron losses and probably 
decrease B, and then the Ib/Ibp would be lower also due to 
the higher temperature of the protons. 

On the other hand, this maximum magnetic field 
appears significantly below equipartition with radiation. 
Our comparison with the data of Cygnus X-1 suggests, 
Ib/Ibh <0.3. Nevertheless it suggests that the emission of 
the corona is not powered by the magnetic field, as assumed 
in most accretion disc corona models (see e.g. di Matteo, 
Celotti & Fabian 1997; Merloni & Fabian 2001). Indeed if 
the corona is powered by the magnetic field, then the lu- 
minosity can be written as L = V^C/s/tdis where tdia is the 
magnetic field dissipation time scale. Since L = !7r V^/iesc, 
we have Ub/Ur = tdis/^csc. Our constraint on the magnetic 
field in Cygnus X-1 then translates into; 

fdis< 0.3(1 + rT/3)7?/c~0.45i?/c (14) 

which points toward extremely fast magnetic dissipation. 
The magnetic dissipation time-scale is expected to be tdis ~ 
R/va with va < c/3 is the Alfven speed. This would gives 
idis ^ 3R/c, in contradiction with equation 1141 

We stress however that this requirement of a low mag- 
netic field in the LHS is entirely due to the observed strength 
of the non-thermal MeV excess in the spectrum of Cygnus 
X-1 and our interpretation of it as the non-thermal tail of 



the Comptonising electron distribution. If contrary to our 
one-zone hypothesis, the MeV excess is produced in a differ- 
ent region than the bulk of the thermal Compton emission, 
the constraint might be relaxed. 

5.2.2 High soft state 

In the HSS the mag netic field (2.3 x 10® G for our assumed 
size of the emitting region) is 6 times larger than in the LHS, 
with a magnetic energy density exceeding both the particle 
and radiation energy density. The magnetic field must be 
higher in the HSS, in order to keep the temperature of the 
thermalised electrons high enough (see Section [4.2.31) . Such 
a magnetic field would be strong enough to power the corona 
and this supports the magnetic corona model for the HSS. 
We note however that spectral fits of other HSS spectra 
of Cygnus X-1 by G99, do not require any heating of the 
thermal electrons {hh = 0), in which case the magnetic field 
would be much lower. Actually, constraining the magnetic 
field in the HSS would require a more detailed comparison 
with a larger set of data. 

5.3 Coulomb heating by hot protons 

Our model including Coulomb heating by hot protons pro- 
vides a good description of the spectra of Cygnus X-1, al- 
though in both states Coulomb heating is not enough and 
some level of non-thermal acceleration is required to explain 
the non-thermal tails. 

5.3.1 Low hard state of Cygnus X-1 

In the LHS, both our simulations and analytical estimates 
suggest a proton temperature which is rather low when com- 
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pared to what obtained in typical two-temperature solu- 
tions. The low proton temperature we infer is related to 
the measured Thomson optical depth and luminosity. The 
optical depth in Cygnus X-1 is larger than unity and there- 
fore the electrons are quite efflciently coupled with the pro- 
tons. Hot protons would produce a luminosity that is larger 
than observed. We estimated an upper limit on the pro- 
ton temperature at about 25 MeV, possibly much lower if 
there are other forms of electron heating or acceleration. 
In contrast, typical ADAF solutions have kTi in the range 
10-100 MeV. Having a temperature comparable to ADAF 
models would require a size of the emission region more than 
5 times smaller than the value (10* cm) assumed in this 
work. In Cygnus X-1, the hard state emission is well mod- 
elled as coming from a hot inner flow surrounded by a cold 
disc truncated at 30-100 -Rg, (e.g., DGK07). If so, having 
a high proton temperature would imply a mass of the black 
hole lower than 4Mq which is clearly excluded. 

5.3.2 Low hard states of XTE J1118+480 and GX 339~4 

We cannot exclude however that a combination of low mass, 
smaller size (in terms of gravitational radii) and uncertain- 
ties related to our assumption of an homogeneous spherical 
region, may combine to underestimate the proton temper- 
ature in Cygnus X-1 by a factor of 5. However, there are 
other sources which have a much lower luminosity and simi- 
lar Comptonisation parameters. For instance Frontera et al. 
(2001b) analyse the BeppoS'^Xspectrum of XTE J1118-I-480 
during its outburst of 2000, while the source had a luminos- 
ity of about one order of magnitude lower than Cygnus X-1. 
The observed X-ray flux of XTE J1118-)-480 was 4.8x10'' 
erg cm~^ s~^, which injected in equation (|13p and @ gives 
Ti/Tc < 1 + 25. Fixing the black hole mass at its minimum 
possible value (6 Mq) and largest distance allowed by the 
uncertainties (d=1.8±0.6 kpc, McClintock et al. 2001), and 
keeping i?=30 Rg, one obtain a conservative upper limit on 
Ti/Ts <10 which is much lower than what is expected in 
two-temperature accretion flows. Then reaching high pro- 
tons temperatures would require a very small R < 3-Rg- 
Similar conclusions could be reached for the source GX 339- 
4 which was observed at comparable flux levels, with similar 
electron temperature and optical depth (see e.g. Joinet et 
al. 2007). For lower luminosity sources, it becomes difficult 
to obtain robust constraints on the optical depth. At least, 
for sources above 10~^ Le a two-temperature plasma seems 
quite unlikely. 

5.3.3 Comparison with standard two -temperature 
accretion flow models 

Yuan & Zdziarski (2004) show that two-temperature ac- 
cretion flow models tend to produce higher electron tem- 
peratures and lower optical depths than what is typically 
measured in Cygnus X-1, GX 339-4 or XTE J1118-H480. 
This is in agreement with what we find here: a high proton 
temperature would indeed require a Thomson optical depth 
lower than what is observed. Also among the simulations 
shown in Yuan and Zdziarski (2004), the models in which 
the plasma parameters are the closest to the observed ones 
have a luminosity much larger than that of GX 339-4 and 
XTE J1118-I-480 in their low luminosity states. 



Esin et al. (1998) presents a comparison of the ADAF 
model spectra with Ginga and CGRO / OSSE data of Cygnus 
X-1 in the LHS. In this comparison, Esin et al. assume a 
black hole mass of M=9Mq and a distance d=2.5 kpc, which 
(for ii=30 Rg) implies 5 =3.5. They found a relatively good 
agreement for the shape of the spectrum. However, although 
they assumed a large 5, their model luminosity is still a 
factor of three larger than what is observed (they had to re- 
normalise their model spectra to fit the observations) . This is 
in agreement with the estimates presented here: with 5 — 3.5 
a luminosity larger by a factor of 3 would allow for a proton 
temperature larger than 10 MeV. 

Esin et al. (2001) show that the simultaneous Ghan- 
dra/RXTE spectra of XTE J1118-f480 is qualitatively 
matched by the ADAF model. Unfortunately the ion tem- 
perature and the plasma parameters are not given in this 
paper and we cannot understand this result. However we 
note that the shape of the model spectra present significant 
differences with the data at high energy, which suggests that 
the temperature and optical depth of the ADAF are different 
from the accurate value measured by Frontera et al. (2001) 
with BeppoSAX. Similarly, Yuan, Cui & Narayan (2005) 
compare the same data with a coupled ADAF/jet model 
(in which the X-ray emission is dominated by the ADAF) 
and find a good agreement. Nonetheless, in their best fit 
model, the region where most of the radiation is produced 
has a Thomson optical depth tt of only 0.2 and the electron 
temperature is around 600 keV (F. Yuan private communi- 
cation). Unlike the BeppoSAX data, the RXTE data do not 
allow one to constrain the electron temperature. Because of 
their lack of sensitivity around 100 keV, the RXTE data 
allow for optically thin high temperature solutions that are 
naturally produced in ADAF models but would probably be 
excluded by the BeppoSAXdaXa.. 

In any case, our estimates for the proton temperature 
also relies on the poorly constrained size of the emitting 
region and we can not completely exclude that this dimen- 
sion is much smaller than 30-Rg. Actually, the emissivity of 
an ADAF and of other two-temperature accretion flows de- 
pends on the distance to the black hole (or radius) . Our size 
of the emitting region R should correspond to the radius 
(or range of radii) where most of the luminosity is radiated. 
Yuan & Zdziardski (2004) show that, at least at high lu- 
minosity, the emission rate per logarithm of radius peaks at 
about 3-4 Rg- This suggests a much more compact emitting 
region than the one considered here. This could allow such 
models to match the spectra and luminosity of Cygnus X-1 
and possibly XTE J1118-f480. An investigation of whether 
and how this is possible, would be very interesting but is 
largely out of the scope of the present work. 

We therefore conclude that in the HSS of black hole bi- 
naries with accurately measured optical depth (tt —1.5-3) 
and electron temperature (fcTo ~50-100 keV), a hot two- 
temperature accretion flow is viable only if the X-ray emis- 
sion is concentrated in a very compact region (lower than a 
few -Rg)- 

5.3.4 High soft state 

In contrast in the HSS of Cygnus X-1 the optical depth is 
lower than in the LHS, and a much higher proton temper- 
ature (~50 MeV) is allowed by the energy balance. This 
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would be in agreement with two-temperature accretion disc 
coronae models that have been proposed as alternative to 
magnetically dominated coronae (di Matteo, Blackman & 
Fabian 1997; Rozahska & Czerny 2000). However, from the 
power-law shape of the high energy HSS spectra most of 
the power must be released to the electrons through non- 
thermal processes. Moreover, as mentioned in Section 15.21 
it is not clear whether thermal heating is really required in 
the HSS and the protons could be cold as well. 

5.4 Slope of the injected electron distribution and 
acceleration processes 

5.4-1 Constraints on the injection slope of the electrons 

In the HSS, the slope of the injected electrons Finj is di- 
rectly constrained by the shape of the high energy pow- 
erlaw spectrum (in the Thomson limit, the photon index 
r = rinj/2 + 1). Our rough comparison allows us to con- 
strain Finj in the range 2.1-2.7 depending on the details 
of the model. In the LHS on the other hand the slope of 
the non-thermal electron distribution is not very well con- 
strained by the shape of the weak high energy excess. How- 
ever, if injection of non-thermal particles is the only form 
of heating for the plasma, then this slope must be much 
steeper (Finj > 3.5) in order to limit the synchrotron emis- 
sion (which otherwise would cool down the thermal compo- 
nent). Even if some form of thermal heating of the electrons 
(through Coulomb heating for instance) can maintain the 
high temperature of the electrons, the injection slope must 
be steep in order to keep the SSC luminosity at the observed 
level. This is a very strong constraint on the injection mech- 
anism in the LHS. 

5.4.2 Consequences for acceleration models 

The different injection slopes in the HSS and LHS could 
be due to different acceleration mechanisms, but they could 
also be due to the same physical acceleration process oper- 
ating in a different regime (e.g. different magnetic and pho- 
ton fields, electron density...). Investigating this issue would 
require a more detailed modelling of the acceleration pro- 
cesses. Our results indicate that neither magnetic dissipation 
nor Coulomb heating in a standard two-temperature plasma, 
is the main heating/acceleration process of the corona in the 
hard state. 

There is also another important issue regarding the ac- 
celeration process: It is not clear whether injecting non- 
thermal electrons with a power-law distribution provides a 
good description of the acceleration processes at work in 
the corona of black hole binaries. For instance, Belmont et 
al. (2008b) show that when the synchrotron boiler effects 
are taken into account, stochastic acceleration processes in 
a magnetised plasma lead to an electron distribution which 
is very close to a Maxwellian. Departure from a Maxwellian 
occurs only in some specific conditions when the stochastic 
acceleration process can accelerate particles only above a 
given threshold energy. A more detailed modelling and inves- 
tigation of these acceleration processes and their relevance 
to the observations is deferred to a future work. However, 
the results of Belmont et al. (2008b), and, in a different con- 
text, those of Katarzyhski et al. (2006), already show that it 



is difficult to avoid a large fraction of the power being trans- 
ferred directly to the thermal electrons. If non-thermal ac- 
celeration processes indeed lead to a direct thermal heating 
of the electrons, this would reduce room for a contribution 
from Coulomb heating in the LHS. For instance, this would 
certainly reduce dramatically the uppers limits on the pro- 
ton temperature that we have set for Cygnus X-1 and pos- 
sibly exclude the presence of a two-temperature plasma in 
this source. On the other hand modelling the high soft state 
with such acceleration models might be more challenging 
and should set serious constraints on the acceleration mech- 
anisms. 



6 CONCLUSION 

We have explored the effects of the synchrotron self- 
absorption on the particle energy distribution and emitted 
radiation of the corona of accreting black holes. We have 
found that, in both spectral states of black hole binaries the 
coronal emission could be powered by a similar non-thermal 
acceleration mechanism. In the LHS the synchrotron and e-e 
Coulomb boilers redistribute the energy of the non-thermal 
particles to form and keep a quasi-thermal electron distri- 
bution at a relatively high temperature, so that most of the 
luminosity is released through quasi thermal comptonisa- 
tion. In the HSS, the soft photon flux from the accretion 
disc becomes very strong and cools down the electrons, re- 
ducing the thermal Compton emissivity. Then most of the 
emission is produced by disc photons up-scattered by the 
non-thermal cooling electrons. Another difference between 
the two states is that the slope of the injected electrons has 
to be steeper in the LHS to reduce the synchrotron emission. 

Our comparison of simulations with the high energy 
spectra of Cygnus X-1 in the LHS allowed us to set up- 
per hmits on the magnetic field and the proton tempera- 
ture. Our results indicates that the magnetic field is at most 
at equipartition with the particle energy density and below 
equipatition with radiation (unlike what is assumed in most 
accretion disc corona models). The proton temperature is 
found to be lower than 25 MeV. 

The constraint on the proton temperature is indepen- 
dent of the synchrotron boiler effects (it is valid even if the 
magnetic field B=0). Indeed, in a two-temperature plasma, 
due to electron ion coupling, there is a relationship between 
the luminosity, the size of the emitting region, the ion and 
electron temperatures, and the Thomson optical depth. For 
the measured luminosity, electron temperature and optical 
depth of several black hole binaries, like XTE J1118-f480 or 
GX 339-4, a proton temperature much higher than the elec- 
tron temperature Ti/Te ~100-1000 (as in a two-temperature 
accretion flow) requires an extremely small size of the emit- 
ting region (a few _Rg or less). 

In the HSS the constraints on the magnetic held and the 
proton temperature are less stringent. The magnetic field 
could power the emission of the corona. 

However these constraints on the plasma parameters 
rely on qualitative comparisons of the simulated spectra 
with the data. It would be interesting to refine and set more 
quantitative constraints on the magnetic field, the proton 
temperature and other physical parameters of the corona 
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of accreting black holes. We hope to present such detailed 
spectral fits in a forthcoming work. 
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APPENDIX A: ANALYTICAL ESTIMATES OF 
THE STEADY-STATE NON-THERMAL 
ELECTRON DISTRIBUTION AND 
SYNCHROTRON TURN-OVER FREQUENCY 



The Cyclotron frequency, 
eB 



(Al) 



can be expressed in terms of the parameters of our model as 
follows: 



Ec — hvc/rUcC = 



3AZb 



(A2) 



where A is the Compton wavelength of the electron and 
is the fine structure constant. 

For a relativistic electron of Lorentz factor 7 most of the 
synchrotron power is emitted around v ~ Therefore, 
for a wide range of electron energy distributions, the syn- 
chrotron emission at a frequency v is dominated by electrons 
of Lorentz factor: 
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(A3) 



where e = hv/m^(? . 

For non-thermal electrons the radiative losses can be 
written as follows: 
4c 



7r 



(A4) 



where g accounts for the Compton losses. It is easy to see 
that g is larger at large 1/Ib- It is however rather difficult to 
estimate accurately. Because most of the radiation energy 
density is in the form of X-ray photons, the Compton inter- 
actions with the non-thermal electron occur mainly in the 
Klein-Nishina regime. As a consequence g is considerably 
reduced with respect to its expected Thomson limit value. 
For the models with inth/^B = 3, the radiation and mag- 
netic field energy density are comparable and yet the losses 
of the non-thermal particles are dominated by synchrotron 
{g ~ 0.1). Another complication is that in the Klein-Nishina 
limit g depends slightly on 7 (i.e. the Compton losses do not 
scale linearly with ■y'^P'^). Since in most of our models, syn- 
chrotron dominate the losses of the non-thermal electrons, 
we can neglect this dependence and a rough estimate of g is 
enough. 

We found that a reasonable approximation of g could 
be obtained by neglecting the losses on the Comptonised X- 
ray photons. If one consider as targets only the synchrotron 
and external soft photons in the Thomson regime: 



3 ^sync ~\~ 



47r 



where 



Lsync^ T 



(A5) 



(A6) 



and I/syiic is the self-absorbed synchrotron luminosity of the 
electrons. An estimate for /sync is derived below which then 
allows to evaluate g (see equations (|A20[I and (|A21|) ). 

On the other hand the e-e Coulomb losses can be esti- 
mated as: 

3c IuAtt 



7c 



2R f3{j) 



(A7) 



where {7) is the average Lorentz factor of the targets elec- 
trons (see Nayakshin & Melia 1998) and is close to unity (we 
use {7)^1.15 in our numerical evaluations). And therefore: 
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At high energies, the losses are dominated by Compton and 
synchrotron while at low energies Coulomb dominates the 
cooling. The transition from ratiative to Coulomb cooling 
occurs at a reduced momentum pc such that: 



Pc = C (Pc + 1), 

where 

^ _ 9 In Att 



8/b(1 + <?){7)' 



(A9) 



(AlO) 



At large ^, pc — \^ while at small ^, pc — For 
the cases, Zb=0.1, 10, and 1000, this gives pc=19.4, 1.98 and 
0.56 respectively. ^ characterises the relative importance of 
e-e coulomb as a cooling process. Coulomb collisions become 
less and less important at larger compactness. 

Using the 'stationary' approximation (see e.g. Fabian 
et al. 1986, Ligthman & Zdziarski 1987, Coppi 1992), the 
steady state non-thermal electron distribution can be esti- 
mated as: 



-(7) = 



d'y f(7)/(7r +7c). 



(All) 



where T{'y) is the distribution of injected electrons. Assum- 
ing 7niin = 1 and Finj > 2, for 1 < 7 < 7max, the non- 
thermal particle distribution can be estimated as 

r(7) = ro7^"^'°-ip~^/(l +C7P"''), 
with 



(A12) 



ro = -|-^(l+p)-^F(rinj), 
167r (B 



and 
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Fin 
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(A13) 



(A14) 



This expression is in good agreement with the results of our 
simulations for 7 > 7t and when the Compton losses are not 
too strong (i.e. Inth/is ^ 10). At lower energies synchrotron 
self-absorption effects become important and the stationary 
approximation gradually breaks down. 

At energies above Pc, the distribution is a powerlaw 
r(7) — Toy'" with s = 1 -I- Finj. For such a distribution the 
reduced synchrotron self-absorption frequency et is given by 
(see WZOl): 



et = 



3{3+s) 



Tg Gl 



l/(4+s) 



J2+s)/(4+s) 



(A15) 



where Gi ~ 1 is defined in equation (9) of WZOl. 

Then the Lorentz factor of the electrons emitting 
around the turn over frequency, 71 = y^^-, can be writ- 
ten as. 



7i = 




(A16) 
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If the particles emitting around the turn-over frequency 
are in the coulomb dominated regime ( i.e. the evaluation 
of equation (|A16|) leads to 74 < y'pf+1) then Coulomb 
effects must be taken into account: 

^(2+s)/(4+.)^ (A17) 

with s now equal to the power-law index in the Coulomb 
regime s = Finj - 




7* 
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(A18) 



According to both equations (|A16|) an (|A18[) . 7t is insensitive 
to the model parameters. In the regime studied in this work 
7t ^ 10. 

Now that we have an estimate of 7t we can evaluate the 
synchrotron compactness (as defined by equation (|A6|) ): 



hync ^ -^Ib I r(7)7^/3^d7. 



(A19) 



In the case where 7t > yjji + l and in the relativis- 
tic limit, the synchrotron radiation compactness, controlling 
the Compton cooling of the Maxwellian component, can be 
estimated as: 



Is 



Inth -F(rinj) 2-ri, 

T7t 



Combining equation (|A20|I with ([ASp . and setting 7t = 10, 
we can evalutate g as: 
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(A21) 



Using this estimate of g, a better approximation of 7t can 
be obtained using equation (|A16p . which in turn may be 
used improve the estimate of g and so on. Since 7t depends 
only weakly on g, convergence is usually achieved after a 
few iterations only. For the general case a similar iterative 
procedure can be applied numerically using equation ()A19p . 

The synchrotron boiler plays a role in the heating of the 
thermalised particles if, at 7t, synchrotron dominates over 
Coulomb cooling: 



2r,3 ^ 91nArT 
8 iBiV 



(A22) 



Since 7t ~ 10, the synchrotron boiler is important for 
thermal heating if Ib/tt 1/5. 



This paper has been typeset from a TjjK/ I^T^^X file prepared 
by the author. 
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